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Abstract We have manufactured more than 250 nom- 
inaUy identical paraffin-coated Cs vapor ceUs (30 mm 
diameter bulbs) for multi-channel atomic magnetome- 
ter applications. We describe our dedicated cell charac- 
terization apparatus. For each cell we have determined 
the intrinsic longitudinal, Iqi, and transverse, I02, relax- 
ation rates. Our best cell shows Ioi/2tt « 0.5 Hz, and 
7o2/27r w 2 Hz. We find a strong correlation of both 
relaxation rates which we explain in terms of reservoir 
and spin exchange relaxation. For each cell we have de- 
termined the optimal combination of rf and laser pow- 
ers which yield the highest sensitivity to magnetic field 
changes. Out of all produced cells, 90% are found to 
have magnetometric sensitivities in the range of 9 to 
30 fT/VHz. Noise analysis shows that the magnetome- 
ters operated with such cells have a sensitivity close to 
the fundamental photon shot noise limit. 



1 Introduction 

Spin polarized alkali vapors prepared by optical pump- 
ing have been used since 50 years for fundamental stud- 
ies in atomic physics and applications thereof [T]. The 
achieved sensitivities depend mainly on the (transverse) 
lifetime T2 of the spin coherence in the vapor, and, to a 
lesser extend, on the (longitudinal) lifetime Ti of the spin 
polarization. Those lifetimes are related to correspond- 
ing relaxation rates by = T^^- To assure a long-lived 
spin polarization, the vapor cells are either filled with 
a buffer gas mixture or are left evacuated while apply- 
ing an anti-relaxation coating on the walls. In the first 
case, the buffer gas in the cell confines the atoms to a 
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diffusion- limited volume and thus reduces the rate of de- 
polarizing wall collisions. In the second case, a thin film 
of paraffin or similar substance applied to the cell wall 
reduces the collisional sticking time with the wall and 
thereby the dephasing interactions with magnetic impu- 
rities embedded in the walls. 

Alkali vapors in paraffin-coated cells were introduced 
in 1958 ^ and have since been widely applied in atomic 
physics spanning applications from magnetometers [31IH 
H], over slow hght studies [6], to spin-squeezing [7], and 
light-induced atomic desorption (LIAD) [51[5] studies. 

Our group develops atomic magnetometers for the 
accurate measurement of small changes in already weak 
fields (typically 10% of the earth's field [TO]), a technique 
that we currently apply to the measurement of the faint 
magnetic fields produced by the beating human heart) 
[lT | [T2 P T3] and for magnetic field measurement and con- 
trol in the search for a neutron electric dipole moment 
[T4lll5j. Both experiments call for a large number (50 to 
100) of individual sensors to be operated simultaneously. 
Although buffer gas cells were used in our initial work 
[TT] , we currently focus on paraffin-coated cells that have 
a reduced sensitivity to magnetic field gradients because 
of motional narrowing and to temperature effects com- 
pared to buffer gas cells P^ll7j . 

In order to fulfill the requirements of the mentioned 
experiments we have initiated a large scale production 
of cells that has yielded over 250 cells in the past year. 
We have developed an automatic cell characterization 
facility for determining the quality and reproducibility 
of the cell coatings. In this work we describe this char- 
acterization facility in detail and report results (intrin- 
sic relaxation times, intrinsic magnetometric sensitivity) 
based on significant cell statistics. A comparative study 
of a small sample of paraffin-coated cells produced over 
four decades was reported in |18j . To our knowledge our 
present study involves the largest sample of coated cells 
ever compared. 
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Fig. 1 (Color online) ParafRn-coated Cs vapor cell. The 
small amount of the solid alkali metal is well visible in the 
sidearm. The arrow points to the capillary which reduces de- 
polarizing collisions of vapor atoms from the cell with the 
solid Cs. 



2 Cell production 



3 X 10^^ mbar, and at the end of filling the cells are 
pumped down to a pressure below 10^^ mbar before be- 
ing sealed. The filled cells are activated by heating them 
in a oven at 80 °C for 10 hours, while assuring that the 
sidearm is kept at a sightly lower temperature. In this 
way we produce 10 coated cells in one week. 

3 The cell characterization setup 

Following manufacture, each cell undergoes a characteri- 
zation procedure in a dedicated experimental apparatus 
for determining the relevant parameters that indicate 
its magnetometric properties. Our current magnetome- 
ters use the technique of optically detected magnetic 
resonance in the Double Resonance Orientation Magne- 
tometer or DROM configuration (notation introduced in 
[12]), also called M^j-configuration [l9lfT4] . The underly- 
ing theory will be addressed below. It was thus a natural 
choice to use the same technique for the dedicated cell 
testing facility. 



The paraffin-coated glass cells are manufactured in our 
institute. Pyrex is formed into a spherical bulb (inner 
diameter of w 28 mm, wall thickness of 1 mm) that is 
connected to a sidearm consisting of a Pyrex tube with 
4 mm inner (7 mm outer) diameter, which acts as a 
reservoir to hold the droplet of solid cesium after coating, 
filling, and sealing the cell (Fig.[T]). The metallic Cs is the 
source for the saturated Cs vapor filling the cell. Near 
the cell proper, the sidearm is constricted into a capillary 
with a design diameter of 0.75(25) mm that reduces spin 
depolarizing collisions with the bulk Cs in the sidearm. 

A typical coating and filling process takes about one 
week. Ten cells are mounted on a glass structure together 
with a paraffin containing reservoir and a Cs metal con- 
taining ampule, both isolated from the vacuum system 
by break-seals. The system is connected to a turbomolec- 
ular pump stand via a liquid nitrogen cold trap and all 
coating and filling steps are performed in a vacuum be- 
low 10~^ mbar. Prior to coating, the whole structure is 
baked for 5 hours at 370 °C. 

The coating process is similar to the one reported in 
[51121 ■ Our current choice of coating material is a com- 
mercial paraffin, Paraflint Hi, from Sasol Wax Ameri- 
can Inc. After baking the system, the break-seal of the 
paraffin reservoir is broken by a piece of iron sealed in 
a glass bead ( "hammer" ) manipulated from the outside 
by a permanent magnet. The wax is deposited onto the 
cell walls by heating the paraffin reservoir. During the 
coating procedure the pressure rises to 9 x 10~^ mbar, 
and the cell is kept isolated from the cesium containing 
ampule. Once the cell is coated, the same hammer is 
used to break the seal of the Cs ampule and a thin film 
of metallic Cs is distilled into the cell's sidearm by heat- 
ing the Cs ampule, after which the end of the sidearm 
is sealed off. During Cs distillation the pressure rises to 



3.1 Experimental setup and signal recording 

The experimental apparatus is shown in Fig. 121 The laser 
source is a DFB laser (A = 894 nm) whose frequency is 
actively stabilized to the 4 => 3 hyperfine component 
of the Cs Di transition using the dichroic atomic vapor 
laser lock (DAVLL) technique [2(1. The light is carried 




Fig. 2 The cell testing apparatus. Frequency stabilized laser 
light is carried by a multimode fiber into a threefold magnetic 
shield (L: lenses). Circular polarization is created by a polar- 
izer (P) and a quarter- wave plate (A/4). The transmitted 
power is recorded by a photodiode (PD) and the modulated 
light power components are extracted by a lock-in ampli- 
fier. A personal computer controls the light power, performs 
scans of the frequency u via a programmable frequency syn- 
thesizer (PFS, Stanford Research model SR345), and records 
the lock-in signals. 
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by a 400 /iin diameter multimode fiber into a three-layer 
mu-metal magnetic shield that contains the actual dou- 
ble resonance setup. Prior to entering the fiber, the laser 
power, Pi, is computer-controlled via a stepper- motor 
driving a half-wave plate located before a linear polar- 
izer. The light leaving the fiber is coUimated and passes 
a linear polarizer followed by a quarter-wave plate to 
create circular polarization before entering the Cs cell. 
The fiber is wound into several loops so that the exit- 
ing light is completely depolarized, thus avoiding vibra- 
tion related polarization fluctuations that translate into 
power fluctuations after the polarizer. 

The paraffin-coated Cs cell to be characterized is 
placed in the center of the magnetic shields where three 
pairs of Helmholtz coils and three pairs of anti-Helmholtz 
coils compensate residual stray magnetic fields and gra- 
dients, respectively. A static magnetic field Bq with an 
amplitude of a few fiT is applied in the yz-plane at 45° 
with respect to the laser beam direction, k = z. The 
transmitted light power is recorded by a nonmagnetic 
photodiode and then amplified. Absorbed laser light pumps 
the Cs atoms into the nonabsorbing (dark) |F=4, Mp—i, 4) 
magnetic sublevels, thereby creating a vector spin polar- 
ization (orientation) (x (Fz). A small magnetic field 
rf-field Bi{t) of a few nT, constant in amplitude, but 
rotating at frequency w, is applied in the plane perpen- 
dicular to Bq. The choice of a rotating, rather than a lin- 
early polarized, oscillating field is used to suppress mag- 
netic resonance transitions in the F—3 state pTI. Bi{t) 
drives magnetic resonance transitions between adjacent 
sublevels in the F—A hyperfine state, whose Zeeman de- 
generacy is lifted by the static magnetic field B^. For 
a properly oriented magnetic field Bq the transmitted 
light power will be modulated at the rotation frequency 

UJ. 

When Lu is close to the Larmor frequency ul — 7_f-Bo, 
where ~ 2n ■ 3.5 Hz/nT is the Cs ground state gyro- 
magnetic factor, a resonance occurs in the absorption 
process, manifesting itself in both the amplitude and 
phase of the light power modulation. The corresponding 
in-phase, quadrature, and phase signals are extracted 
by means of a lock-in amplifier (LIA, Stanford Instru- 
ments, model SR830) whose output signals are read by 
a personal computer. The rotating field frequency is gen- 
erated by a computer controlled programmable synthe- 
sizer. The computer varies this frequency, w, by a linear 
ramp in the range of ± 27r • 100 Hz around the Lar- 
mor frequency during a scan time of 40 s. A dedicated 
electronics box generates from this AC voltage two 90° 
dephased AC currents that drive two perpendicular coil 
pairs (not shown in Fig. [2|) producing the rotating field 
Bi{t). 

The characterization of each individual cell consists 
in the recording of resonance spectra for a set of 12 
selected (and computer controlled) values of the laser 
power Pl in the range of 1 to 12 /iW. It is difficult to 
determine the absolute laser intensity for a given laser 



power Pi, because of the (asymmetric) transverse beam 
profiles and their modification by the cell's spherical 
shape. We therefore quantify the light intensity in terms 
of the laser power P^ , to which it is proportional. Note 
that Pl used below refers to the power measured af- 
ter the cell with the laser frequency resonant with the 
4-^3 Cs Di transition and the rf power off. A typical au- 
tomated characterization run, including insertion of the 
cell into the apparatus, takes 10 minutes. Data analysis is 
performed by a semi-automatic dedicated Mathematica [22] 
code, which takes another 5 minutes. In a regular work- 
ing day it is thus possible to characterize 30 to 40 cells. 



3.2 DROM theory 

A modulation of the transmitted power only occurs when 
the static magnetic field Bq is neither parallel nor per- 
pendicular to the direction of light propagation. In that 
case the transmitted light power has components that 
oscillate in phase, Dui^ and in quadrature, A^j, with re- 
spect to the rotating field 



B,{t)^^e 
If 



(1) 



The in-phase and the quadrature components depend 
on the detuning, 5 = w — wq, between the driving, w, and 
the Larmor, wq, frequencies. The dependence of and 

on i5 are dispersive and absorptive Lorentzians given 

by m 



D^{S) = -77 (P^) sin (20) 
AUS)^-V{F.) sin {29) 



(2) 



where Aq = ?/ {Fz) sin (26) is a common signal amplitude 
that depends — via the spin polarization {Fz) created by 
optical pumping and the detection of the polarization's 
precession via light absorption — on the laser power P^ . 
The calibration constant 77 includes all of the apparatus 
constants such that D^{5) and At^[5) are measured in 
Vohs. With respect to Fig. H U{t) = D^{5)cosujt + 
A^{S) sin ujt. The phase (pujiS) between the drive and the 
power modulation 



I2 

J ((5) = + arctan ( — 



(3) 



depends also on the detuning S. The expressions ^ and 
([3]) are valid for atomic media with an arbitrary ground 
state angular momentum, as may be shown easily by 
a theoretical treatment analogous to the discussion of 
the signals in the DRAM (double resonance alignment 
magnetometer) geometry presented in |12j . In the above 
expressions, 6 is the angle between the applied magnetic 
field Bq and the laser beam propagation direction k. 
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Fig. 3 (color online) Lock-in demodulated magnetic reso- 
nance signals. Top: The dispersive signal (blue) represents 
the in-phase component D{u}) and the absorptive signal (red) 
the quadrature component A{ijj). Bottom: Phase signal 4>{'-^)- 
Experimental points are shown together with lines fitted ac- 
cording to ((ll-Q. All signals were recorded at Bo — 4 f/T 
(tjo ^ 27r ■ 11640 Hz), with Pz, = 6 ^iW, and Bi = 1.3 nT. 

3.3 Signal analysis 

Since the resonance signals are extracted by a lock-in 
amplifier, and since it is experimentally difficult to pre- 
cisely determine the phase of the rotating field (and 
hence the phase difference between that field and the 
modulation of the photocurrent), the signals produced 
by the lock-in amplifier are superpositions of the absorp- 
tive and dispersive lineshapes A^[5) and D^{5). Using 
the fitting procedure described in detail in [lU] it is pos- 
sible to extract the pure absorptive and dispersive com- 
ponents. For fitting the theoretical lineshapes the com- 
bined apparatus constants Aq = rj {Fz) sin {29) is taken 
as one fitting parameter, with Aq measured in Volts. 
Other parameters are the relaxation rates !{ and I2 , the 
resonance frequency ljq, an unknown overall phase, as 
well as weighting factors of the absorptive and disper- 
sive components. The Rabi frequency I2rf can be easily 
calibrated as described in 15] and a fixed numerical value 
is used when fitting ^ and 

Typical resonance lineshapes of the in-phase, quadra- 
ture, and phase signals are shown in Fig. [31 together with 
the fitted theoretical shapes ([2]) and ^ . Fitting the ab- 
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Fig. 4 Laser power dependence of the relaxation rates li 
(boxes) and I2 (diamonds). The experimental points are fit- 
ted with (|4|. The (statistical) error bars on the individual 
data points are smaller than the symbol size. 

sorptive and dispersive spectra by ^ with the relax- 
ation rates 11 and I2 as free parameters yields a strong 
correlation between the two rates in the x^-minimizing 
algorithm, with corresponding large uncertainties in the 
numerical values. We have therefore opted for the follow- 
ing fitting procedure. In a first step, we use the fact that 
the phase does not depend on 7i and fit the dependence 
0(u;) given by ([3]) to the data. The resulting L2, value is 
then used as a fixed parameter in the subsequent simul- 
taneous fit of the absorptive and dispersive lineshapes 
to infer 11. In this way we obtain (/f, 72)-pairs for each 
value of the laser power P^. In addition, the fits yield 
the overall signal amplitude Aq. 

4 Results 

^.1 Relaxation rates 

Figure H] shows the dependence of the longitudinal and 
transverse relaxation rates on the laser power Pl . There 
is, to our knowledge, no theoretical algebraic expression 
describing that dependence for ground states of arbitrary 
angular momentum F. We therefore fit, as in the 
dependence by a quadratic polynomial 

^Pl) ^o^ + a, Pl + 13^ Pi , (4) 

which allows us to infer the intrinsic relaxation rates, Tqi 
and 7o2, i.e., the relaxation rates extrapolated to zero 
light power. 

4.2 Signal amplitudes 

Figure [5] shows the dependence of the signal amplitude 
^0 on the laser power P^. Here again, we have no theo- 
retically derived algebraic expression describing that de- 
pendence for transitions between states with arbitrary 
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angular momenta. We therefore, as in [3], fit the experi- 
mental dependence by the empirical saturation formula 



So{Pl)=C 



Pi 



{Pl + Psi){Pl + PS2) 



(5) 



which accounts for an amplitude growing as P| at low 
powers, and where Pgi and Ps2 are saturation powers. 

Figures, m and [5] show typical dependencies of ii, F2, 
and 5*0 on P^ for a given cell, together with the fits (solid 
lines) by ([S]). We have characterized 253 parafhn-coated 
cells of equal diameter using the method described above. 
The histograms in Fig. [6] (top, middle) show the distri- 
butions of the intrinsic longitudinal and transverse re- 
laxation rates of the 241 best cells. The scatter plot in 
the lower graph of Fig. [5] shows that the two rates are 
strongly correlated. The fitted line represents a linear 
relation of the form -/02 = s ^ii + » with s ~ 1.00(1) and 
a = 1.35(3) Hz. The longitudinal and transverse relax- 
ation rates are thus equal, up to a constant offset that 
affects the 7o2 values only. For an isotropic relaxation 
process, in which all Zeeman sublevels relax at the same 
rate, one would expect Tqi — Io2- In section [5] below we 
will come back to a quantitative discussion of those con- 
tributions. 



4.3 Magnetometric sensitivity 

The intrinsic relaxation rates are well suited to char- 
acterize each individual cell. In particular, the trans- 
verse rate _Zo2, which determines the intrinsic width of 
the signals and D^^, is relevant for magnetometric 
applications. However, the intrinsic rates are, by defini- 
tion, rates for vanishing laser and rf powers. Therefore, 
the magnetometric sensitivity of a given cell can not be 
inferred directly from the intrinsic rates, since magne- 
tometers have to be operated at finite laser and rf power 
levels. 
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Fig. 5 Magnetic resonance amplitude versus laser power. 
The experimental points are fitted with the polynomial ex- 
pression from ([5j, which yields, for this specific cell, the sat- 
uration parameters Psi ~ 634 nW and Ps2 = 16.3 ^W. The 
error bars are smaller than the plotting symbol size. 




Foi / 2n (Hz) 

Fig. 6 Histogram of intrinsic longitudinal (top) and trans- 
verse (middle) relaxation times of 241 coated cells. The upper 
axis in the top graph gives the radius of the effective depo- 
larization spot that models reservoir relaxation (see text). 
The lower graph shows the correlation between the relax- 
ation rates, together with a fit of the form 7o2 ~ sloi + a. 



The linear zero crossing of the dispersive signal D^^ 
near resonance is convenient for magnetometric applica- 
tions since any magnetic field change 6B yields a signal 
change 



dD.. 



dB 



5B 



(6) 
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that is proportional to SB. The lowest magnetic field 
change SB that can be detected depends on the shot 
noise of the DC photocurrent cx Pl. A feedback resis- 
tor, Rp, in the transimpedance amplifier, marked I /U in 
Fig. [21 transforms that photocurrent into a photovoltage 
Ul, whose shot noise (in a bandwidth of 1 Hz) is given 

by 



SUl 



Rf SIl = Rf \/2elL Rf 



2QEPLe^ 



hu 



(7) 



where Qe = 70% is the quantum efficiency of the pho- 
todiode, and u the laser frequency. The experimentally 
measured signal noise lies « 20% above the shot noise 
level, due to laser power fluctuations and amplifier noise. 
With the calibration constant 77 in ([2|) , 5D^ is expressed 
in Volts, i.e., in the same units as 5Ul- 

For each set of the experimental parameters Pl and 
Q-cf one can thus define the magnetometric sensitivity as 
the field fluctuation (5i?NEM that induces a signal change 
5Di^ of equal magnitude than 5Ul ■ This noise equivalent 
magnetic field fluctuation (NEM) is thus given by 



6B 



SUr 



NEM 



I dB 



p2 



7F 



^0 



SUl 



f 



(8) 
(9) 



Ao, ii, and I2 are {Pl dependent) parameters obtained 
from the fits of the experimental D^^ spectra. SU l is as- 
sumed to be the Pl dependent shot noise value ([7]). We 
recall that fi^{ is not a fit parameter, and that calibrated 
numerical values of f^i-f are inserted in ^ when evalu- 
ating ^Bnem- 

For each cell we have evaluated SBj^em for a range of 
parameters Pl and J7rf- Figure [7] shows a typical result 
in terms of a contour plot of 5_Bnem- For each cell we 
determine the optimal NEM value, ^^neM' by a numer- 
ical minimization procedure. The minimum for the cell 
shown in Fig. [7] is indicated by a cross. 

The distribution of minimal NEM values, SB'^^^, 
thus obtained is represented in form of a histogram in 
Fig. E Only cells with (5-B™m < 4GfT/Vlfe are shown. 
This set represents 94% of all cells we have produced to 
date. 



5 Discussion 

The distribution of linewidths shown in Fig. [6] reveals 
a dependence of the form Iq2 = Iqi + Z\77ciax, with a 
constant offset relaxation rate Z\i7oiax, whose numeri- 
cal value (fit parameter a in Fig. [6]) is Z\i7ciax/27r — 
1.35 Hz. Here we show that -Zqi is ultimately limited 
by atoms escaping to the sidearm, and that Z\7i^eiax is 
mainly determined by spin exchange collisions {Apix) 
with a minor contribution from magnetic field inhomo- 
geneities {AF^b)- 
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Fig. 7 Plot of (5-Bnem as a function of the amplitude Orf of 
the rotating field and of the laser power Pl- The contours 
represent the lines of constant NEM, spaced by 1 fT/VHz, 
with selected numerical values indicated. The cross refers to 
the minimal NEM value, which, for the cell represented here 
has a value of 10.5 fT/Vllz- 
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Fig. 8 Histogram of the minimal NEM values, ^Bne'mj of 
241 cells, which represent 94% of the cells produced to date. 



5.1 Longitudinal relaxation 

The intrinsic longitudinal relaxation rate Iqi is limited 
by processes which thermalize the magnetic sublevel pop- 
ulations, such as atoms escaping through the capillary 
to the sidearm where they eventually collide with the 
solid Cs droplet, atoms hitting an imperfectly coated 
surface spot of the spherical bulb, or atoms being ab- 
sorbed by the coating 8,9J. AH of those processes can 
be parametrized in terms of an effective depolarizing 



surface area Udc 



' dop- 



We will refer to such pro- 



cesses in general as "reservoir losses" . The distribution 
of Tqi values in the top graph of Fig. [5] represents the 
statistical distribution of such imperfections, due to un- 
controlled parameters in the cell production process. In 
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a spherical cell of radius R the rate of wall collisions is 
7waii = 3v/AR, where v is the average thermal velocity. 
The intrinsic longitudinal relaxation rate can thus be ex- 
pressed in terms of the effective depolarizing spot radius, 
r-dcp, via 



is given by 



3 V rl^p 



iBl=7wall^-^^. (10) 

The upper axis in the top graph of Fig. [5] shows the 
radius rjop corresponding to the Iqi value on the lower 
axis. The best cell produced so far has a longitudinal 
relaxation rate Ioi/27t « 0.50(5) Hz, which corresponds 
to c?dop = 2rdop = 1 mm. This value is compatible with 
the design diameter, dcap = 0.75(25) mm, of the capil- 
lary, which shows that Iqi is ultimately limited by atoms 
escaping into the sidearm. 



5. 2 Transverse relaxation: field inhomogeneities 

If the offset magnetic field Bq varies over the cell volume 
it produces a distribution of resonance frequencies col, 
and hence a broadening of the magnetic resonance lines 
given by ([2]) and fS]). The fitting analysis interprets this 
broadening as an increase of the transverse linewidth Iq2 
by an amount AF^b ■ A main advantage of coated cells 
over buffer gas filled cells is that, because of multiple 
wall collisions, the atoms explore a large fraction of the 
cell volume during the spin coherence time, which effec- 
tively averages out field gradients. Standard line narrow- 
ing theory '23] predicts that an inhomogencous magnetic 
field gives a lowest order contribution 



(11) 



to the transverse relaxation rate, where Z\i?rms is the 
rms value of the magnetic field averaged over the cell 
volume, and Tc the correlation time of the field ffuctu- 
ations seen by the cell, which can be approximated by 
the mean time between wall collisions. This expression 
is valid in the so-called good averaging regime [23j, i.e., 
for 'ypABj-ynsTc ^ 1. From the geometry of the used coils 
we estimate Z\i3rms to be on the order of 2 nT, which 
yields Av^b = FabI'^t^ = 30 mHz. Even when allowing 
for a 5 times larger inhomogeneity (i.e., AB = 10 nT) 
from uncompensated residual fields — recall that we ac- 
tively compensate linear field gradients — one still has 
AvAB < 0.1 Hz. We can thus ascertain that the contri- 
bution from field inhomogeneities to Iq2 is negligible. We 
note that the good averaging conditions for AB = 2 nT 
and 10 nT read 7i?Z\i?r,ns''c = 0.004 and 0.02, respec- 
tively. 



5.3 Transverse relaxation: spin exchange 

As derived by Ressler et al., [21], the contribution from 
spin exchange collisions to the transverse relaxation rate 



An^ = ^ 



21 



21+1 



nCsCTc^Vr , 



(12) 



where / is the nuclear spin, rtcs the Cs number density, 
Vr the relative velocity of colliding atoms, and CTox = 
2.06 X 10"^'' cm^ [24j the spin exchange cross section for 
Cs-Cs collisions. The parameter a describes the slowing 
down of the spin relaxation due to the hyperfine interac- 
tion. In small magnetic fields a « 0.63 for the M= — 4 
M= - 3 transition (Fig. 3 of [24]). At T = 20(1) °C the 
contribution of spin exchange collisions to io2 evaluates 
to 

AT 

-^ = 1.6(2) Hz, (13) 

where the error reflects the uncertainty in the number 
density. This value is compatible with the experimen- 
tal value Iq2 — lii — (27r)1.35(3) Hz. We are therefore 
confident that dephasing spin exchange collisions give 
the main contribution to the transverse relaxation rate, 
notwithstanding a certain scatter of the spin exchange 
cross sections in the literature. 



5. 4 Fundamental limits of magnetometric sensitivity 

The ultimate sensitivity of the type of magnetometers 
described here is limited by two fundamental processes, 
viz., photon shot noise limit and spin projection noise. 
One can show (Appendix [X| that the minimal NEM im- 
posed by the shot noise of the detected photons is given 

by 



rnPSN 
"^NEM 



1 



1 



hv 



2^/2i2 



(14) 



where Pl is the power detected after the cell, Qe the 
quantum efficiency of the photodiode for photons of en- 
ergy hv, and kq the resonant absorption coefficient of 
the driven hyperfine component for unpolarized atoms. 
For a time interval of t ^ 0.5 s, the result corresponds 
to a measurement bandwidth of 1 Hz. In {F^) is the 
spin polarization 



A/=-4 



(15) 



in the i^=4 state, where the P4,m are the populations of 
the magnetic sublevels |F=4, M). The analyzing power 
for the transition F-^F' , Aff', depends in general on 
the applied laser power and accounts for population ef- 
fects such as hyperfine pumping. Its value has been de- 
termined by a numerical model based on rate equations 
[25] . It is a slowly varying function in the domain of laser 
powers considered here, with value A43 = 1.15(5). 
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For our apparatus, kqL « 0.7, Qe — 0.7, so the above 
can be rewritten as 



ArPSN 

"-°NEM 



(fT) 



0.146 



(16) 



For our best cell, A43 (F^) = 0.39(4), ii/27r = 3.40 Hz, 
and I-z/^TT = 4.75 Hz at the optimum laser power ol 
3.6 /LtW, which yields an expected sensitivity of (5i?NEM — 
7.0(7) fT, to be compared with the measured minimal 
NEM of the ceU of 9(1) fT. For a more typical cell with 
Ii/2TT = 10 Hz, ri/2TT = 8.65 Hz, and A43 (F^) = 0.46(5) 
at the optimal power of 5 /iW, the expected minimal 
NEM is SB^%^ = 9.6(1.0) IT, indicating that the shot 
noise limited NEM grows less than linearly in I^. 

Spin projection noise limits the magnetometric sen- 
sitivity to 

^ ' (17) 



rnSPN 

"-Onem 



If 



where A'at = ^PatKeii is the number of atoms in the 
F=4: state that contribute to the signal, with pat be- 
ing the total Cs number density, and T4g11 the cell vol- 
ume. For a measurement time tmeas of 0.5 s, one finds at 
T = 20(1) °C, SB^^li = 0.74(2) fT for i2/27r = 4.75 Hz. 
In our magnetometers spin projection noise thus has a 
negligible contribution. 



maximizing (F^) via hyperfine repumping, which could 
win, at most, a factor of 2-3. 

It is well known that in the spin exchange limited 
regime an increase of the atomic density by heating the 
cell does not increase the magnetometric sensitivity, since 
both I2 and kq in (|14p grow proportionally to the density. 
The same holds for I2 and A^atom in (HZ]) - However, when 
operating the magnetometer in a regime where spin ex- 
change is not the limiting factor, one expects an improve- 
ment of the sensitivity by increasing the atomic number 
density. 

We will use the cells in multi-sensor applications in 
fundamental and applied fields of research. Since an op- 
timal magnetometric sensitivity is reached with a typical 
light power of approximately 5 fiW, a single diode laser 
can drive hundreds of individual sensors 13J. This scal- 
ability, together with the very good reproducibility of 
the coated cell quality reported here, will allow us to re- 
alize in the near future a three-dimensional array of 25 
individual sensors for imaging the magnetic field of the 
beating human heart, a signal with a peak amplitude 
100 pT [13] . With a reliable and inexpensive multichan- 
nel heart measurement system, magnetocardiograms can 
be measured in a few minutes, times which are of interest 
in the real world of clinical applications. 



6 Summary and conclusion 

We have manufactured and characterized a set of 253 
paraffin-coated Cs vapor cells of identical geometry (15 mm 
radius spheres), 90% of which have an intrinsic trans- 
verse relaxation rate in the range of 2 to 6 Hz. Under 
optimized conditions of laser and rf power those cells 
have intrinsic magnetometric sensitivities, 6B 



NEM' - 



I the 

range of 9 to 30 fT/\/Hz under the assumption of (light) 
shot-noise limited operation in a DROM-type magne- 
tometer. 

The magnetometric sensitivity is determined by the 
intrinsic transverse relaxation rate, which, for the best 
cell of our batch has a value of 27r • 2 Hz, of which 
« 0.5 Hz are due to reservoir (Ti) relaxation, and 
« 1.5 Hz are due to spin exchange relaxation. Improv- 
ing the relaxation properties by reducing reservoir re- 
laxation is technologically demanding, and would only 
marginally improve the overall sensitivity. Spin exchange 
relaxation, on the other hand, cannot be suppressed in 
coated cells, although it was shown that spin exchange 
relaxation can be suppressed in high pressure buffer gas 
cells, yielding sub-fT magnetometric sensitivity [26]. We 
thus conclude that our cells are as good as coated cells 
of that diameter can be, disregarding a possible 25% re- 
duction of 7o2/27r by a suppression of reservoir losses. 

The expected photon shot noise limited NEM of our 
cells is very close to the measured NEM. The most promis- 
ing improvement in sensitivity is expected to come from 



A Photon shot noise Hmit 

Consider a light beam of power Pin traversing a vapor 
of thickness L. The transmitted power detected by a 
photodiodc with quantum efficiency Qe is given by 



Pdetit) = QeP„ 



,-K{t)L _ 



QePl 



(18) 



The time dependent absorption coefficient K,{t) consid- 
ering the in-phase component of the magnetic-resonance 
induced modulation is 



n{t) = Ko 1 - AfF' (F,) 



(52 



71 ■^^rf 



cos Ujt 



(19) 

where kq is the resonant optical absorption coefficient 
for a sample of unpolarized atoms and the polarization 
(Fz) is as defined by (fTS)) . The analyzing power for tran- 
sition F-^F' , ApF', accounts for population effects such 
as hyperfine pumping: more details are given in the dis- 
cussion in the main text following (|15p. Lock-in detection 
extracts from (fTS|) the rms value 



pLIA 
^det 



V2 



QE/'^„e-««^ k^LAff' {F, 



(20) 

of the in-phase component of the power modulation. 

Light power P can be converted to a photon count 
N seen during time t via N — for v the photon fre- 
quency. Hence N^J/^ represents the number of photons 
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carrying magnetometric information, thus 



dB 

which evaluates to 



dNg{^ dPt{^ dS 
dPti^ dS dB' 



(21) 



dN 



LIA 

det 



1 



dB 




assuming an rf amplitude {f2,.{ = y/Iill) which 
mizes the result. 

The photon shot noise limited magnetometric sensi- 
tivity is given by 



coPSN 
"^NEM 



dNt{^ 
dB 



where 



Assembling the above components gives 
SB 



-KqL 



PSN ^ 2V2i^ /IT 1 1 / 

IF y riKoLAFF'{F,)\ QEPLt' 



(23) 



(24) 



(25) 



which is the required result. 
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